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Abstract. We calculate the accretion disc temperature profiles, disc luminosities and boundary layer luminosities 
for rapidly rotating neutron stars considering the full effect of general relativity. We compare the theoretical values 
of these quantities with their values inferred from EXOSAT data for four low mass X-ray binary sources: XB 1820- 
30, GX 17+2, GX 9+1 and GX 349+2 and constrain the values of several properties of these sources. According 
to our calculations, the neutron stars in GX 9+1 and GX 349+2 are rapidly rotating and stiffer equations of state 
are unfavoured. 
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1. Introduction 

A low mass X-ray binary (LMXB) is believed to contain 
either a weakly magnetised neutron star or a black hole as 
the central accretor. The X-ray emission arises from the 
innermost region of the accretion disc around the compact 
star. In the case of a neutron star, there is an additional 
X-ray component coming from the boundary layer of the 
star. Mitsuda et al. (1984) showed that the spectrum of 
a luminous LMXB can be fitted by the sum of a single 
temperature blackbody spectrum (believed to come from 
the boundary layer) and a multicolour blackbody spec¬ 
trum (may be originated from the accretion disc). However 
these authors used Newtonian models to fit the observed 
spectra. But near the surface of a neutron star, the ac¬ 
cretion flow is expected to be governed by the laws of 
general relativity due to the presence of strong gravity. 
Therefore general relativistic models should be used for 
the purpose of fitting to get the correct best-fit values of 
the parameters. Besides, the principal motivation behind 
the study of the spectral and temporal behaviours of neu¬ 
tron star LMXBs is to understand the properties of very 
high (~ 10 15 g cm -3 ) density matter at the neutron star 
core (van der Klis 2000). This is a fundamental problem 
of physics, which can not be addressed by any kind of lab¬ 
oratory experiment. The only way to answer this question 
is to assume an equation of state (EOS) model for the neu¬ 
tron star core, to calculate the structure parameters of the 
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neutron star and hence to calculate an appropriate spec¬ 
tral model. By fitting such models (for different chosen 
EOSs) to the observed data, one can hope to constrain the 
existing EOS models and hence to understand the prop¬ 
erties of high density matter. However, general relativistic 
calculation is essential to calculate the structure parame¬ 
ters of a neutron star and therefore to constrain the EOS 
models. 

It is expected that the neutron stars in LMXBs are 
rapidly rotating due to accretion-induced angular momen¬ 
tum transfer. LMXBs are thought to be the progenitors of 
milli-second (ms) radio pulsars (Bhattacharya & van den 
Heuvel 1991) like PSR 1937+21 with P ~ 1.56 ms (Backer 
et al. 1982). The recent discovery of ms (P ~ 2.49 ms) X- 
ray pulsations in XTE J1808-369 (Wijnands & van der 
Klis 1998) has strengthened this hypothesis. Therefore it 
is necessary to calculate the structure of a rotating neutron 
star considering the full effect of general relativity. This 
was done by Cook, Shapiro & Teukolsky (1994) and the 
same procedure was used by Thampan & Datta (1998), 
to calculate the luminosities of the disc and the boundary 
layer. 

The disc temperature profile for a rapidly rotating 
neutron star was first calculated by Bhattacharyya et al. 
(2000). These authors also compared their theoretical re¬ 
sults with the EXOSAT data (analysed by White, Stella 
& Parrnar 1988) to constrain different properties of the 
LMXB source Cygnus X-2. The present work is a continu¬ 
ation of theirs, in which we constrain several properties of 
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four LMXB sources: XB 1820-30, GX 17+2, GX 9+1 and 
GX 349+2, using the same procedure. These sources were 
also observed by EXOSAT and the data were analysed by 
White et al. (1988). 

XB 1820-30 is an atoll source which shows type I X- 
ray bursts. GX 17+2 and GX 349+2 are Z sources, of 
which the former shows X-ray bursts. GX 9+1 is an atoll 
source. As all of them are LMXBs (van Paradijs 1995), 
the magnetic field of the neutron stars are believed to be 
decayed to lower values (~ 10 8 G; see Bhattacharya & 
Datta 1996 and Bhattacharya & van den Heuvel 1991). 
Therefore, we ignore the effect of the magnetic field on 
the accretion disc structure in our calculations. 

In section 2, we give the formalism of the work. We 
present the results and discussion in section 3 and give a 
summary in section 4. 


2. Formalism 

2.1. Theoretical formulae 

In order to calculate the disc temperature profile, the disc 
luminosity and the boundary layer luminosity for a rapidly 
rotating neutron star considering the full effect of general 
relativity, we need to compute the structure of the ro¬ 
tating star. To do this, following Cook et al. (1994), we 
choose a stationary, axisymmetric, asymptotically flat and 
reflection-symmetric (about the equatorial plane) metric, 
given by 

dS 2 = g^dx^dx” (fi,v = 0,1,2,3) 

= -e 7+p dt 2 + e 2a (df 2 + f 2 d0 2 ) 

+e 7_p f 2 sin 2 Q(d</> — ujdt) 2 , (1) 

where the metric coefficients 7 , p, a and the angular speed 
(co) of zero-angular-momentum-observer (ZAMO) with re¬ 
spect to infinity, are all functions of the quasi-isotropic ra¬ 
dial coordinate (f) and polar angle (9). The quantity r is 
related to the Schwarzschild-like radial coordinate (r) by 
the equation r = fe^ 7_p ^ 2 . We use the geometric units 
c = G = 1 in this paper. 

With the assumption that the star is rigidly rotating 
and a perfect fluid, we solve Einstein’s field equations and 
the equation of hydrostatic equilibrium (the equations are 
given in the appendix) self-consistently and numerically 
from the centre of the star upto infinity to obtain the 
metric coefficients and II* (angular speed of neutron star 
with respect to infinity) as functions of r and 9. The inputs 
of this calculation are a chosen EOS and assumed values 
of the central density and the ratio of polar to equatorial 
radii. The outputs are bulk structure parameters, such 
as gravitational mass (M), equatorial radius ( R ), angular 
momentum (J), moment of inertia (I) etc. of the neutron 
star. We can also calculate the specific disc luminosity 
(Ed), the specific boundary layer luminosity (-Ebl), the 
radius (r or b) of the innermost stable circular orbit (ISCO) 
and specific energy ( E ), specific angular momentum ( l) 
& angular speed (Hk) of a test particle in a Keplerian 


orbit (see Thampan & Datta 1998 for a description of the 
method of calculation). 

Then we calculate the effective temperature of a geo¬ 
metrically thin blackbody disc, which is given by 

T eS (r) = ( F(r)/aY / 4 (2) 


where <7 is the Stefan-Boltzmann constant and F is the 
X-ray energy flux (due to viscous dissipation) per unit 
surface area. We calculate F using the expression of Page 
& Thorne (1974): 


. M , . 

F (r) = ~ A — d(r) 

47rr 

where 

q(r) = -n K ,r(E - 


(E — fl^l)l r dr 


( 3 ) 

( 4 ) 


Here ri n is the disc inner edge radius and a comma fol¬ 
lowed by a variable as subscript to a quantity, represents 
a derivative of the quantity with respect to the variable. 
The values for E and l are given by the two conditions (cir¬ 
cularity and extremum) for orbits (see Thampan & Datta 
1998; Bhattacharyya et al. 2000): 


E — u)l = 
l = 


e (7+P)/2 
a / 1 — V 2 

vfe ^~ p ^ 2 

a / 1 — V 2 


( 5 ) 

( 6 ) 


where v = (f2 — co)fe~ p sin9 is the physical velocity of the 
matter. The equation of motion in the azimuthal direction 
and that in the time direction yield the Keplerian angular 
speed as 


D k 


c 2 p(f) ~ l /r 2 

(E - col) 


u(r) 


( 7 ) 


Eq. (4), that is exactly valid for a black hole accre- 
tor, is also valid for a neutron star accretor, if the disc 
does not touch the star. But, for a disc that extends upto 
the surface of the neutron star, the torque at the disc in¬ 
ner edge will not vanish, and eq. (4) will not be strictly 
valid for such a case. However, for a very rapidly rotating 
neutron star, the angular speed of a particle at the disc 
inner edge will be close to the spin rate of the star, and 
hence the torque is expected to be negligible. According 
to Bhattacharyya et al. (2000), even when the spin rate 
is not large, eq. (4) should give approximately correct re¬ 
sults. This is because, Page & Thorne (1974) argued that 
the error in the calculation of T e g will not be substan¬ 
tial outside a radial distance l.ln n . In our calculation, we 
find that almost whole of the disc X-ray flux comes from 
well-outside this radial distance. 

Eq. (2j) gives the effective disc temperature T e g(r) with 
respect to an observer comoving with the disc. For our pur¬ 
pose, this expression of the temperature must be changed 
to that in the observer’s frame, taking into account the 
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gravitational redshift and the rotational Doppler effect. In 
order to keep our analysis tractable, we use the expression 
given in Hanawa (1989) for this modification : 


With this correction for (1 + z), we define a temperature 
relevant for observations {T^ x ) as: 

T imax ^ rrimax /ro 

obs ^ 7— J eft 
t + Z 

where the superscript ‘max’ denotes the maximum tem¬ 
perature value in the profile. T)£ ax is related to the colour 
temperature (T’™f x ) by T^ x = T c ™ 1 ax //, where f is 
the colour factor (Shimura & Takahara 1995; see also 
Bhattacharyya et al. 2000 for details). We compare the 
observationally inferred value of the maximum tempera¬ 
ture with T™f x . It is to be remembered that eq. (8) is 
valid for a non-rotating neutron star and a ‘face-on’ (i.e., 
inclination angle i = 0) observation. However, the error in¬ 
volved by these assumptions is expected to be well-within 
the error bars considered in subsection 2.3. 

The structure of a neutron star for a given EOS is 
described uniquely by two parameters : the gravitational 
mass (M) and the angular speed (fl*). For each adopted 
EOS, we construct constant M equilibrium sequences with 
fl* varying from the non-rotating case (static limit) upto 
the centrifugal mass shed limit (rotation rate at which in¬ 
wardly directed gravitational forces are balanced by out¬ 
wardly directed centrifugal forces). So we are able to cal¬ 
culate T e I g ax , Ed and Ebb as functions of M and Q* for a 
chosen EOS model. 

2.2. Equations of State 

The neutron star structure parameters are quite sensitive 
to the chosen EOS models. For the purpose of a general 
study, we have used the same four EOSs as considered in 
Bhattacharyya et al. (2000), namely, (A) Pandharipande 
(1971), (B) Baldo, Bombaci & Burgio (1997), (C) Walecka 
(1974) and (D) Sahu, Basu & Datta (1993). EOS model 
(A) is for hyperonic matter. It is assumed that hyper- 
onic potentials are similar to the nucleon-nucleon poten¬ 
tials, but altered suitably to represent the different isospin 
states. Model (B) is a microscopic EOS for asymmet¬ 
ric nuclear matter, derived from the Brueckner Bethe- 
Goldstone many- body theory with explicit three-body 
terms. The three-body force parameters are adjusted to 
give a reasonable saturation point for nuclear matter. EOS 
model (C) corresponds to pure neutron matter and is 
based on a mean-field theory with exchange of scaler and 
(isoscalar) vector mesons representing the nuclear inter¬ 
action. Model (D) is a field theoretical EOS for neutron- 
rich matter in beta equilibrium based on the chiral sigma 
model. The model includes an isoscalar vector field gener¬ 
ated dynamically and reproduces the empirical values of 
the nuclear matter saturation density and binding energy 
and also the isospin symmetry coefficient for asymmetric 


nuclear matter. Of these, model (A) is soft, (B) is inter¬ 
mediate in stiffness and (C) & (D) are stiff EOSs, with 
(D) as the stiffest. 

2.3. Constraining procedure 

We choose four LMXB sources observed by EXOSAT 
(data analysed by White et al. 1988). For each source, we 
take the best-fit values of the parameters T™ ax , Ld and 
Lbl, where L denotes the luminosity. On the other hand, 
for the chosen values of /, the accretion rate (M) and 
M, we theoretically calculate the values of T™f x , Ld and 
Lbl as functions of O*, for an assumed EOS model. Then 
comparing these theoretical values with the observed ones, 
we constrain different parameters of the chosen source 
(see Bhattacharyya et al. 2000 for detailed description). 
However, to take into account the uncertainties in the 
fitting procedure and in the value of z, and also those 
arising due to the simplicity of the model, we consider a 
range of acceptable values for T“ 1 ax , Ld and Lbl- We take 
two combinations of deviations around the best-fit values, 
namely, (10%, 25%) and (20%, 50%), where the first num¬ 
ber in parentheses corresponds to the error in T c “ ax and 
the second to the error in the best-fit luminosities. 

3. Results and Discussion 

In this paper, we calculate gravitational mass sequences 
for different EOS models and constrain several properties 
of four LMXB sources. For the neutron star in each of 
the sources, we assume M =1.4 M 0 (i.e., the canonical 
mass value for neutron stars). We take two values for cosi 
(i is the inclination angle of the source) for each source, 
namely, 0.2 and 0.8. These two widely different values en¬ 
sure the sufficient generality of our results. For the four 
sources, the best-fit values (White et al. 1988) of the pa¬ 
rameters T™f x , Ld and Lbl are given in Table 1. 

We take the distance ( D ) of the source XB 1820-30 as 
6.4 kpc (Bloser et al. 2000). We assume D = 8 kpc for both 
GX 17+2 and GX 9+1, as their locations are believed to 
be near the galactic centre (Deutscli et al. 1999; Hertz et 
al. 1990) and distance of the galactic centre is 7.9 ± 0.3 
kpc, as concluded by McNamara et al. (2000). For GX 
349+2, we take D = 9 kpc (Deutsch et al. 1999). 

We display the constrained values with the help of four 
tables. It is to be noted that here M is presented in unit 
of M e = 1.4 x 10 17 M/Mq g s -1 . The Eddington rate is 
M e /q, with q = Ebl + Ed- Therefore, as the actual value 
of q is much lesser than 1.0 (generally not greater than 0.3 
and for rapidly rotating neutron star, typically less than 
0.2), the value of Eddington accretion rate is much higher 
than M e . For all the sources, as the stiffness of the EOS 
models increases, the absolute values of the allowed spin 
frequencies (i/*) and rotational frequencies in the ISCO 
(z^in) decreases. This is because, for a stiffer EOS model, 
neutron star radius is higher and it can support lesser 
amount of rotation. The energy conversion efficiency is 
also lesser for a stiffer EOS model (as the neutron star 
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for this case is less compact) and therefore higher accre¬ 
tion rate is needed to generate the observed luminosity 
(as seen from the tables). In the following, we describe 
the results for four sources in four subsections and give a 
general discussion in the last subsection. 

3.1. XB 1820-30 

We display the allowed ranges of different parameters for 
the source XB 1820-30 in Table 2. It is seen that for cos * = 
0.2, the spin frequency (y+j of the neutron star comes out 
to be very high. But in the case of cost = 0.8, it is not 
possible to constrain (*/*) for (20%,50%) uncertainty set 
(for all EOS models) and for both the uncertainty sets 
(for EOS model D). The ranges of the colour factor are in 
general consistent with the results of Shimura & Takahara 
(/ ~ 1.7 — 2.0) and Borozdin et al. (1999) (/ = 2.6). 
However, some discrepancy can be noted with the latter 
one for softer EOS models & cos * = 0.2. The value of the 
rotational frequency in the ISCO (*q n ) comes out to be ~ 1 
kHz for all the cases. The values of the stellar equatorial 
radius are in the range 8 — 21 km. The peak of the disk 
effective temperature occurs in the radial range 18 — 30 
km, and always well-outside (by several kilometers) the 
neutron star’s surface. This shows the validity of eq. (4), 
as discussed in section 2. The overall range of the accretion 
rate comes out to be 0.5 — 31.4M e . 

3.2. GX 17+2 

Table 3 shows the results for the source GX 17+2. Here the 
ranges of (*/*) are similar to those for XB 1820-30. But for 
GX 17+2, the value of * is expected to be moderately high 
(Titarchuk et al. 2001) and therefore cos* is not possibly 
as high as 0.8. It is, therefore, quite likely that the neutron 
star in this source is rapidly rotating, although no decisive 
statement can be made. The allowed values for / for GX 
17+2 is systematically lower than those for XB 1820-30, 
and in the case of softer EOS models & cos* = 0.2 they 
do not tally with the result of Borozdin et al. (1999). The 
allowed values of (zq n ) corns out to be ~ 1 kHz, but for 
EOS model (A) & cos* = 0.8, 2 kHz value is also possible. 
The ranges of R and r^f x are similar to those for XB 
1820-30 and the allowed values of the accretion rate are 
in the range 2.0 — 131.0M e . 

3.3. GX 9+1 

The results for the source GX 9+1 are given in Table 4. 
For this source, *q,-value comes out to be very high for all 
the EOS models and for both the cos*-values. Here, the 
allowed values for / are inconsistent with Borozdin et al. 
(1999) for softer EOS models & cos* = 0.2. The allowed 
values of *q n are ~ 1 kHz and the allowed ranges of R and 
r”fj ax are 10 — 21 km and 18 — 30 km respectively. The 
allowed values of accretion rate for this source come out 
to be in the range 3.8 — 116.7 M e . 


3.4. GX 349+2 

The allowed ranges of different parameters for the source 
GX 349+2 are given in Table 5. As is the case for GX 9+1, 
here also the value of comes out to be very high for 
all the chosen EOS models and cos*-values. The allowed 
values for f in general tally with the results of Shimura 
& Takahara (1995), but like other three sources do not 
match with the result of Borozdin et al. (1999) for softer 
EOS models & cos* = 0.2. Here comes out to be ~ 1 
kHz, and R & r^g x are in the ranges 10 — 21 km and 
18 — 30 km respectively (like GX 9+1). The accretion rate 
for this source is in the range 4.5 — 168.5 M e . 

3.5. General Discussion 

Here we have constrained the values of several properties 
of four LMXB sources. For all of them, the accretion rates 
come out to be very high (always > 0.5 M e ). This is in 
accord with the fact that these are very luminous sources. 

The rotation rate of neutron star in each of the sources 
is very high (close to the mass shed limit) for cos* = 0.2. 
This is because, the values of Tbl/^d are very low for 
these cases (see Thampan & Datta 1998, Bhattacharyya 
et al. 2000). But, for cos* = 0.8, rotation rate can not 
be constrained effectively for the sources XB 1820-30 and 
GX 17+2. Therefore, for these two sources, no general con¬ 
clusion (about the values of H*) can be drawn. However, 
the allowed ranges (combined for all the cases considered 
in a table) of are 0.93 — 1.00 and 0.75 — 1.00 

for the other two sources GX 9+1 and GX 349+2 respec¬ 
tively (here fl^ms is the fl* at the mass shed limit; see 
Bhattacharyya et al. 2000 for the mass shed limit values). 
Therefore the neutron stars in these two sources can be 
concluded to be rapidly rotating in general. 

Our calculated allowed ranges for f are in accord 
with the results obtained by Shimura & Takahara (1995). 
However, if we take the value / = 2.6 (obtained by 
Borozdin et al. 1999), one would require a very stiff EOS 
model or a mass greater than M = 1.4 M@ for most of 
the cases with cos* = 0.2. 

High frequency quasi- periodic-oscillations (kHz QPO) 
have been observed for three (XB 1820-30, GX 17+2 and 
GX 349+2) of the chosen sources. The observed maxi¬ 
mum kHz QPO frequencies are 1.100 kHz (XB 1820-30), 
1.080 kHz (GX 17+2) and 1.020 kHz (GX 349+2) (van der 
Klis 2000). Now, as pointed out in Bhattacharyya et al. 
(2000), the maximum possible frequency (i.e., the shortest 
time scale) of the system should be given by the rotational 
frequency in ISCO (*q n ; col. 5 of the tables). Therefore, 
the stiffest EOS model D is unfavoured for cos* = 0.2 
for the source XB 1820-30, as the maximum value of v\ n 
(= 0.941 kHz, Table 2) is less than the observed maximum 
kHz QPO frequency. For the same reason, EOS model D 
is unfavoured for cos* = 0.2 for the source GX 17+2. It 
can also be seen from Table 3 that if we use only the 
narrower limits on the luminosities and colour tempera¬ 
ture, EOS model D (for cos* = 0.8) and EOS model C 
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(for cos* = 0.2) are unfavoured for the same source. Same 
is true for EOS model D for the source GX 349+2. As 
we also see from Table 5, EOS model C is unfavoured 
for cos?' = 0.2 for this source. Therefore, we may conclude 
that the stiffer EOS models are unfavoured by our results. 

We have ignored the magnetic fields of the neutron 
stars in our calculations. Therefore, the necessary condi¬ 
tion for the validity of our results is that the Alfven radius 
(r\) be less than the radius of the inner edge of the disc. 
This condition will always be valid if R > ?’a holds. Here 
ta is given by (Shapiro & Teukolsky 1983), 


-2/7 

o M . 4/7/ M , *3/7 

r A = 2.9 x 10 8 (-=-) /4o( 10 ) 

M e IVIq 

where M is the mass of the neutron star, /z 30 is the mag¬ 
netic moment in units of 10 30 G cm 3 and r\ is in cm. With 
typical values of the parameters for the chosen sources 
(R = 10 km, M = 1.4Mq and M — 10M e ), the upper 
limit of the neutron star surface magnetic field comes out 
to be about 2 x 10 8 G. Therefore, our results are in general 
valid for the neutron star magnetic field upto of the order 
of 10 8 G. This is a reasonable value for the magnetic field 
of neutron stars in LMXBs, as mentioned in section 1. 

It is also to be noted that our results are valid for a 
thin blackbody disc. However, as the spectra of the sources 
were well-fitted by a multicolour blackbody (plus a black- 
body, presumably coming from the boundary layer; White 
et al. 1988), the assumption of thin blackbody disc may 
be correct. 

4. Summary 

In this paper, we have constrained the values of two neu¬ 
tron star parameters (spin frequency and equatorial ra¬ 
dius) for the four LMXBs: XB 1820-30, GX 17+2, GX 
9+1 and GX 349+2. We have also calculated the allowed 
ranges of the colour factor (for accretion disc), rotational 
frequency {u- m ) of a particle in the ISCO, the radius (r^ ax ) 
where the effective disc temperature is maximum and the 
accretion rate for these sources. These have been done 
for a chosen mass 1.4 M 0 (canonical mass) of the neutron 
stars and two values of inclination angles (cos? = 0 . 2 , 0 . 8 ). 
The whole work has been repeated for four EOS models 
(from very soft to very stiff). 

We have drawn the following main conclusions from 
our study. A comparison between the kHz QPO frequen¬ 
cies (observed from the sources) and our calculated val¬ 
ues of v ln has shown that the stiffer EOS models are un¬ 
favoured. By the constraining procedure, we have got very 
high accretion rates for all the sources, which is in ac¬ 
cordance with their high luminosities. The neutron stars 
in the sources GX 9+1 and GX 349+2 have been found 
to be very rapidly rotating, and those in the other two 
sources may also be rapidly rotating (although we can 
not say decisively). This is in accordance with the belief 
that LMXBs are the progenitors of millisecond pulsars. 


It also shows that while calculating the spectral models, 
taking the rapid rotation of neutron stars into account is 
very important. 

It is difficult to constrain EOS models effectively with 
the poor quality EXOSAT data. However, the present gen¬ 
eration X-ray satellites Chandra and XMM have much 
better resolving power. For example, the resolving power 
of Chandra HETGS is 60 — 1000 in the energy range 
0.5 — 10.0 keV. The future generation X-ray satellite 
Constellation-X will have even better resolving power 
(upto 3000). With the spectral data of these X-ray ob¬ 
servatories, it may be possible to constrain EOS models 
and other parameters of LMXBs effectively. Therefore we 
propose that it is essential to compute EOS dependent 
general relativistic spectral models to utilise these good 
quality data in a fruitful way. 
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Appendix 

Here we give the Einstein’s field equations and the equation 
of hydrostatic equilibrium, that were solved for the computa¬ 
tion of the structure of the rapidly rotating neutron star. For 
an axisymmetric and equatorial plane symmetric configura¬ 
tion, the computational domain in spherical polar coordinates 
covers 0 < r < oo and 0 < 9 < n/2. For numerical convenience, 
we make a change of variables (r —> s and 0 —> p) given by 
f = and 6 = cos - 1 /r, where f e is the quasi-isotropic 

radial coordinate of the equator. It is easy to see that s and p 
vary in the range 0 <s<l& 0 <p<l and at the equator 
s = 0.5. 

The four Einstein’s equations (Cook et al. 1994) to solve 
are given below: 


p(s,M) = 


- _ c -t/2 


^P2n(/i) 


1 — S 


2n+l 


ds's'- 


I o (! - S') 2n+2 


xj dfj,'P 2n {n')Sp(s',n')+yYZ^ 


/- 1 ds'(l - S ') 2n -i f 1 , , n , 

X / -^27Ti- / p 2 n(p )Sp{s ,p. 


( 11 ) 


7 (s,n) = ~ 


2 e 7 ^ 2 sin[( 2 n — 1 ) 0 ] 


1—5 

7 r ( 2 n — l)sin# LV s 

n =1 


ds's ' 2 "- 1 f 1 


I o (l-s ') 2n+1 


dp' sin[( 2 n — l) 0 ']S 7 (s', p') 


+ 


1 — s 


rl ds'(l-s ') 2 


x / dp! sin[(2n — 1 ) 0 / ]S 7 ( 5 / , \x) 
Jo 


( 12 ) 


w(s,p) = 


00 nl / \ 

_ _^{2p- 1 )/2 fJn-VM 
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n =1 
/ l2n 


2 n( 2 n — 1 ) sin 0 
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2n+l 
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Q (1 — s') 2n + 2 

2n — 2 


f dp' sin 6' Pl n _ 1 {n')Sa{s , p' 
do 


1 — s 


ds'(l - s ') 2 


x / dp' sin 9'P 2 n - 1 (n')Sa(s', p') 
do 


(13) 


^(p,m + — {(1 — M )[1 + s (1-s) 7, s ] L 


+ [-M+(1-/07,*.]}" 


-{s(l - s)[s(l - s) 7 ,s],s 


+s 2 (1-s) 2 7 2 s -[(1-M 2 )7,J,m 

_ 7,*.[—/* + (1 - M 2 )7,m]}[-M + (1 - M 2 )7,m] 

+ ^[s 2 (! - s) 2 (p, s + 7,7 - (1 - M 2 )(p,m + 7,7] 
x [-M + (1 - 7)7,*.] - «(1 - s)(l - 7) 

x (^(/>,s + 7, «)(/>,*. +7,*.) +7,s*. +7,*7,**) [1 + S (1 ~ s)7, s ] 
+s(l - s)/* 7 , s [l + s(l — s) 7 , s ] + 1(1 — p 2 )e~ 2p 

x ] 2—^—(1 - /x 2 )d>, s d>, M [l + s(l - s) 7 , 4 ] 

( \ 

+ (1-7)7,*.]}] (14) 

where P„(p) are the Legendre polynomials, P™(p) are the as¬ 
sociated Legendre polynomials and sin(n0) is a function of p 
through 6 = cos -1 p. The effective sources S 's are defined as 
(Cook et al. 1994) 
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where u> = r e ui and fi* = r e fl*. Here the tilde over a vari¬ 
able represents the corresponding dimensionless quantity and 
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the variables e and P (not to be confused with Legendre poly¬ 
nomials, that always have a function of ‘n’ as the subscript) 
represent mass-energy density and pressure respectively. 

The equation of hydrostatic equilibrium for a barytropic 
fluid is 


h(P) - h p 


dP 
i+P 

In u* — In «p — / 

J ,c 



(18) 


where h(P) is the dimensionless specific enthalpy as a function 
of pressure and P p ,u p and h p are the dimensionless values of 
pressure, t-component of the four-velocity and the specific en¬ 
thalpy at the pole. The quantity is the (dimensionless) 
central value of the angular speed, which on the rotation axis 
is constant and equal to its value at the pole. The quantity 
F(£l*) = is obtained from an integrability condition on 

the equation of hydrostatic equilibrium. Choosing the form of 
this function fixes the rotation law for the matter. Following 
Komatsu et al. (1989), we set it to F(Cl t ) = A 2 (C f* jC — f1*). 
Here A is a rotation constant such that rigid rotation is 
achieved in the limit A —> oo. An appropriately chosen value 
of h p defines the surface of the star. 
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Table 1 . Best-fit values of the parameters (see text) for four LMXBs. 


Source name 

mmax 

col 

(10 7 K) 

(10 38 
cos i = 0.2 

Lr> 

ergs s -1 ) 
cos i — 0.8 

£bl 

(10 38 ergs s” 1 ) 

XB 1820-30 

1.59 

1.49 

0.37 

0.26 

GX 17+2 

1.76 

6.49 

1.62 

0.71 

GX 9+1 

2.25 

6.01 

1.50 

0.25 

GX 349+2 

2.07 

8.54 

2.14 

0.48 


Table 2. Observational constraints for various EOS models : (A), (B), (C), (D) for the source XB 1820-30. L and U 
stand for lower and upper limits. The parameters are / (colour factor), iz* (rotational frequency of the neutron star), 
t'in (rotational frequency in the ISCO), R (equatorial radius of the neutron star), r^ ax (radius where the effective 
temperature of the disc is maximum) and M (the accretion rate). The limits are for 25% uncertainty in luminosities 
and 10% uncertainty in the colour temperature. Values in brackets are for 50% uncertainty in luminosities and 20% 
uncertainty in the colour temperature, i is the inclination angle of the source with respect to the observer. The mass 
of the neutron star is assumed to be 1.4 Mq. The accretion rate is given in unit of M e = 1.4 x 10 17 M/Mq g s -1 , 
where M is the neutron star mass. 


EOS cos i f v, v in R M 

kHz kHz km km M e 

JA) 02 L 1.31[1.10] 1.751(1.726] 1.756(1.756] 11.2(10.2] 18.6(18.1] 7.7[3.8] 

U 1.91(2.41] 1.755(1.755] 1.819(2.078] 11.4(11.4] 18.7(18.7] 16.1(20.3] 

(B) 0.2 L 1.45(1.30] 1.103(1.059] 1.137(1.126] 15.0(13.7] 23.0(21.6] 9.7[4.5] 

U 2.07(2.61] 1.112(1.113] 1.197(1.372] 15.5(15.6] 23.5(23.6] 19.4(24.4] 

(C) 0.2 L 1.49(1.30] 0.961(0.913] 0.979(0.973] 16.5(15.0] 24.8(23.1] 10.4(4.9] 

U 2.12(2.71] 0.968(0.968] 1.042(1.206] 17.2(17.2] 25.5(25.6] 21.2(26.7] 

■(D) 02 L 1.59(1.40] 0.735(0.687] 0.748(0.743] 19.9(17.7] 29.1(26.5] 12.2(5.6] 

U 2.25(2.81] 0.740(0.740] 0.795(0.941] 20.6(20.7] 30.0(30.1] 25.0(31.4] 

(A) 0.8 L 1.79(1.50] 1.463(0.000] 1.822(1.571] 9.9[7.5] 18.1(18.1] 1.2[0.5] 

U 3.06(4.70] 1.751(1.754] 2.165(2.165] 11.2(11.4] 20.4(22.3] 4.5[6.4] 

(B) 0.8 L 1.94(1.71] 0.498(0.000] 1.207(1.152] 11.3(11.0] 20.2(20.2] 1.4[0.9] 

U 3.22(4.20] 1.102(1.110] 1.782(1.782] 14.9(15.4] 22.9(23.4] 5.6(7.9] 

■(C) 08 L 1.99(1.70] 0.175(0.000] 1.046(0.991] 12.3(12.3] 21.5(21.0] 1.5[1.0] 

U 3.36(4.02] 0.960(0.966] 1.573(1.573] 16.5(17.0] 24.7(25.4] 6.1[8.5] 

(D) 0.8 L 2.11(1.80] 0.000(0.000] 0.806(0.758] 14.7(14.7] 23.1(23.1] 1.8[1.2] 

_U 3.30(3.90] 0.733(0.739] 1.212(1.212] 19.7(20.5] 28.9(29.8] 7.2[9.9] 


Table 3. Observational constraints for various EOS models : (A), (B), (C), (D) for the source GX 17+2. Other 
specifications are same as in Table 2. 


EOS cos i f v, z/ in R rST M 

kHz kHz km km M e 

JA) 02 L 1.01(1.00] 1.754(1.748] 1.756(1.756] 11.4(11.0] 18.7(18.5] 36.1(19.8] 

U 1.43(1.82] 1.755(1.755] 1.773(1.869] 11.4(11.4] 18.7(18.7] 68.7(82.6] 

(B) 0.2 L 1.12(1.00] 1.108(1.097] 1.128(1.122] 15.3(14.7] 23.3(22.6] 44.4(24.4] 

U 1.57(1.96] 1.113(1.113] 1.163(1.236] 15.6(15.7] 23.6(23.6] 82.6(101.7] 

(C) 0.2 L 1.15(1.00] 0.966(0.954] 0.974(0.971] 16.9(16.0] 25.2(24.2] 47.5(26.1] 

U 1.62(2.01] 0.968(0.968] 1.000(1.091] 17.2(17.2] 25.5(25.6] 88.5(111.5] 

(D) 0.2 L 1.22(1.03] 0.738(0.728] 0.744(0.742] 20.3(19.2] 29.7(28.3] 55.9(30.7] 

U 1.72(2.13] 0.740(0.740] 0.766(0.849] 20.7(20.7] 30.1(30.1] 106.4(131.0] 

(A) 0.8 L 1.39(1.20] 1.702(0.000] 1.782(1.571] 9.9[7.5] 18.1(18.1] 6.6[2.0] 

U 2.20(3.72] 1.754(1.755] 2.166(2.166] 11.3(11.4] 18.6(22.3] 18.5(25.5] 

(B) 0.8 L 1.53(1.31] 1.009(0.000] 1.172(1.141] 13.1(11.0] 21.1(20.2] 7.7[3.2] 

U 2.31(3.35] 1.107(1.111] 1.463(1.782] 15.2(15.5] 23.2(23.5] 22.8(30.7] 

(C) 0.8 L 1.61(1.30] 0.858(0.000] 1.010(0.983] 14.3(12.3] 22.4(21.0] 8.3[3.6] 

U 2.33(3.20] 0.965(0.967] 1.289(1.568] 16.8(17.1] 25.1(25.4] 25.0(32.9] 

(D) 0.8 L 1.66(1.41] 0.631(0.000] 0.775(0.750] 16.8(14.7] 25.4(23.1] 9.5[4.4] 

_U 2.47(3.10] 0.737(0.740] 1.011(1.212] 20.2(20.6] 29.5(30.0] 29.3(38.6] 
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Table 4. Observational constraints for various EOS models : (A), (B), (C), (D) for the source GX 9+1. Other 
specifications are same as in Table 2. 


EOS cosi / v, z/ in R r3f x M 

kHz kHz km km M e 

■(A) 02 L 1.33(1.13] 1.755(1.755] 1.756(1.756] 11.4(11.4] 18.7(18.7] 36.9(22.8] 

U 1.85(2.25] 1.755(1.755] 1.756(1.761] 11.4(11.4] 18.7(18.7] 59.9(72.0] 

(B) 0.2 L 1.47(1.24] 1.112(1.110] 1.120(1.117] 15.6(15.4] 23.6(23.4] 43.4(27.4] 

U 2.04(2.49] 1.114(1.114] 1.130(1.147] 15.7(15.7] 23.6(23.7] 73.6(90.6] 

(C) 0.2 L 1.51(1.28] 0.968(0.967] 0.970(0.970] 17.2(17.1] 25.5(25.4] 47.5(29.3] 

U 2.09(2.56] 0.968(0.968] 0.975(0.986] 17.3(17.3] 25.6(25.6] 80.7(99.3] 

(D) 0.2 L 1.61(1.36] 0.740(0.739] 0.742(0.741] 20.7(20.5] 30.1(29.9] 55.9(34.4] 

U 2.24(2.74] 0.740(0.740] 0.745(0.754] 20.7(20.7] 30.1(30.1] 94.9(116.7] 

(A) 0.8 L 1.84(1.61] 1.752(1.728] 1.756(1.756] 11.2(10.3] 18.6(18.1] 7.9[3.8] 

U 2.69(3.52] 1.755(1.755] 1.815(2.064] 11.4(11.4] 18.7(18.7] 16.5(20.3] 

(B) 0.8 L 2.05(1.80] 1.103(1.064] 1.136(1.126] 15.0(13.8] 22.9(21.7] 9.7[4.6] 

U 2.92(3.72] 1.112(1.113] 1.200(1.361] 15.5(15.6] 23.5(23.6] 19.4(24.4] 

■(C) 0+ L 2.10(1.80] 0.961(0.919] 0.978(0.972] 16.5(15.0] 24.8(23.2] 10.4(5.0] 

U 3.00(3.81] 0.968(0.968] 1.041(1.195] 17.2(17.2] 25.5(25.6] 21.2(26.7] 

(D) 0.8 L 2.24(1.92] 0.734(0.692] 0.748(0.743] 19.8(17.8] 29.0(26.7] 12.2(5.7] 

_U 3.19(4.00] 0.740(0.740] 0.802(0.932] 20.6(20.7] 30.0(30.1] 25.0(31.4] 


Table 5. Observational constraints for various EOS models : (A), (B), (C), (D) for the source GX 349+2. Other 
specifications are same as in Table 2. 


EOS cosi / z/, z/ in R rS ax M 

kHz kHz km km M e 

JA) 02 L 1.12(1.00] 1.755(1.754] 1.756(1.756] 11.4(11.4] 18.7(18.7] 50.9(30.7] 

U 1.55(1.92] 1.755(1.755] 1.756(1.771] 11.4(11.4] 18.7(18.7] 86.4(103.9] 

(B) 0.2 L 1.24(1.04] 1.112(1.109] 1.122(1.119] 15.5(15.3] 23.5(23.3] 61.2(37.7] 

U 1.72(2.11] 1.113(1.114] 1.137(1.159] 15.7(15.7] 23.6(23.7] 106.3(130.8] 

(C) 0.2 L 1.27(1.08] 0.968(0.966] 0.971(0.970] 17.2(17.0] 25.5(25.3] 67.1(40.4] 

U 1.76(2.17] 0.968(0.968] 0.979(1.000] 17.2(17.3] 25.6(25.6] 116.6(140.1] 

(D) 0.2 L 1.35(1.14] 0.740(0.738] 0.742(0.741] 20.6(20.4] 30.0(29.7] 78.8(47.5] 

U 1.88(2.31] 0.740(0.740] 0.748(0.765] 20.7(20.7] 30.1(30.1] 136.9(168.5] 

(A) 0.8 L 1.55(1.34] 1.748(1.678] 1.760(1.756] 11.0(9.7] 18.4(18.1] 10.6(4.5] 

U 2.29(3.10] 1.755(1.755] 1.879(2.147] 11.4(11.4] 18.7(18.7] 23.3(30.0] 

(B) 0.8 L 1.71(1.50] 1.096(0.955] 1.148(1.129] 14.6(12.7] 22.6(20.7] 13.1(5.2] 

U 2.47(3.22] 1.110(1.112] 1.244(1.532] 15.4(15.6] 23.4(23.6] 28.0(36.0] 

1C) 0+ L 1.76(1.53] 0.953(0.798] 0.984(0.974] 16.0(13.8] 24.2(21.9] 14.0(5.6] 

U 2.54(3.30] 0.967(0.968] 1.093(1.350] 17.1(17.2] 25.4(25.5] 30.7(38.6] 

(D) 0.8 L 1.87(1.61] 0.727(0.557] 0.753(0.744] 19.1(16.2] 28.2(24.6] 16.5(6.3] 

_U 2.69(3.40] 0.739(0.740] 0.845(1.069] 20.6(20.7] 29.9(30.1] 36.0(46.4] 




